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Abstract 


A novel method of monotone recognition based on the partitioning of the grid 
into discrete structures isomorphic to binary cubes (called “cube-split” 
technique) was proposed in our recent work, and a theoretical level description 
of two algorithms /algorithmic schemes/ solving this problem was also 
introduced. This paper provides implementation details of those algorithms, as 
well as focuses on the recognition of monotone binary functions with a small 
number of units. 
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1. Introduction 


Let 2/* ,, denote the n-th Cartesian degree of the set 2,4, = {0,1,---,m}: 
ered = {(a1, ag An) |Qj € Em+il =1,-- Nn}, 

or, in other words, £7), ,, is the set of vertices of the n-dimensional (m + 1)-valued discrete grid. 
The total number of wertiCes of £7.41, is equal to (m +1)”. We consider a component-wise 
partial order “<” on £7,,, defined in the following way: for arbitrary vertices a = (a4,°*+, An) 
and b = (b,,:::, by) of £7.44, a precedes b (a < b) if and only if a; < b; for i = 1,---,n. Then, 
(Sh41<) is a partially ordered set; we will use its Hasse diagram for geometrical 
interpretations. 

Ff (%4, X2)°%* Xp): Eh, > {0,1} is called a binary function defined on £7 ,,. We say that 
f (X14, X2,°**, Xn) iS a monotone eoener if for any two vertices a,b of 2? ,,, a= b implies: 
f(a) = f(b). The vertices of 2% ,,, where f takes the value “1”, are called units of the function. 
The set of units of f is usually denoted by Ny. The vertices of 5;,4,, where f takes the value “0” 
are called zeros of the function. a’ € £ Sms is called a lower unit of f, if f(a’) = 1 and f(b) = 
0 for every bE EP ,, less than at. a® € Ey, is called an upper zero of f, if f(a°) = 0 and 
f(b) = 1 for every b € £7 ,, greater than a®. 
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Fig. 1 demonstrates the Hasse diagram of 22, and a monotone binary function f defined 
on 32. Highlighted vertices (4,4,4), (4,4,3), (43.4), (3.4.4), (4.4.2), (4.3.3), (3.4.3), (3.3.4), 
(2,4,4), (4,3,2), (3,3,3), (2,4,3), (1.4.4), (1,4,3) are units of the function, where (4,3,2), (3,3,3), 
and (1,4,3) are its lower units. The rest of vertices of 23 are zeros of the function, and (4,4,1), 
(4,2,4), (3,4,2),(2,3,4), and (0,4,4) are its upper zeros. 

We consider the problem of query-based algorithmic identification/recognition of 
monotone binary functions defined on £/1,,. This problem is initially investigated by V. 
Korobkov and V. Alekseev, but also, consecutively, by many other authors [1-4]. For m = 1, 
this is the case of ordinary monotone Boolean functions defined on the binary cube E”, E” = 
{(@4,°**,@,)|@; © {0,1},i = 1,---,n}. Hansel’s chain-splitting technique of E” [5] is a well- 
known effective tool for monotone Boolean function recognition. The outline of the algorithm is 
as follows: the set of vertices of the binary cube is partitioned into disjoint chains of different 
lengths (there are a total of clr 7I chains in the n-dimensional cube). A key property of the 
Hansel’s chains is that once the function values are known for all the vertices in all the chains of 
length k, then the function values, inferable by monotonicity, are unknown for at most two 
vertices in each chain of the next length k + 2. The maximum number of queries to recognize 


the monotone Boolean function defined on the n-dimensional cube is cn ali cln/ a, 
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Fig. 1. Monotone function defined on 23. 


An extension of this technique to the case of multi-valued grids and monotone binary 
functions is obtained in [2-3]. In [2], V. Alekseev developed the algorithm Up, for recognition of 
a monotone binary function defined on 2x, x,..-k, = 5k, X Fk. XX Fes (Fe, = {01,5 — 
1},i = 1,-:-,n), which, in some sense, tries to generalize G.Hansel’s algorithm. [2] proved that: 





Te. ot 
Tey < flloga(k - DI. 
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Where T(Uo) denotes the complexity of Up, and T(Ugp¢) is the complexity of the optimal 
algorithm Ugy:, k = maxk;. It is also found that: 


T(Uopt) 2 |M| + |N| and T(Ug) S |M| + [logok] - INI, 
where M and N are the 2 sets of vertices in the middle layer area of 24. x,..-k,: 
zs 1 
M= {(@y,+, dn) E Skykg-Kky * 41 +--+a, = [; (kj = 1)]}. 
3 1 
N= {Qu dn) © Beegegnky § Ga to + ay = [DE (ki — 1)| +1}. 
In case of £7), ,, the sets in the two middle layers are: 


" mn 
My = {(@y,, dn) Es: a,tet+a,= [= }, 


No = {ay an) ead 2 Ay te tay oe +1}, 


and consequently, the estimate of complexity of the algorithm Up on £7), , will be: 
T(Ug) S |Mol + Llogzm| - |Nol. 


A recent novel method of the monotone recognition based on a partitioning of the grid 
into discrete structures isomorphic to binary cubes (called “cube-split” technique) is proposed in 
[6], and two algorithms /algorithmic schemes/ solving this problem are also introduced. This 
paper provides implementation details of these algorithms, as well as focuses on the recognition 
of monotone binary functions with a small number of units. 

The paper is organized as follows: Section 2 introduces the cube-splitting technique. 
Section 3 provides the implementation framework of these algorithmic schemes. Section 4 
addresses some particular cases with a small number of units that comes from applications. 


2. The Cube-Splitting Technique 


In this section we introduce the cube-splitting technique of recognition of monotone binary 
functions [6]. Two homogeneous areas inside the £7, ,, are defined in the following way: 
- upper homogeneous area 4, - this is the set of all “upper” elements of 5” ,,, i.e., 
elements with all-coordinate values > m/2; 
- lower homogeneous area H, — this is the set of all “lower” elements, i.e., elements with 
all-coordinate values < m/2. It is clear that: 


m+1 


(7) for odd m, 


Al =|A| = 
| | | | (| + 1)” for even m 


The following results were introduced in [6]. 


(1) =7,, can be split into |A | disjoint discrete structures isomorphic to binary cubes: 
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mei = €1 UU Eig, 


where every €; contains exactly one vertex from A, while the remaining vertices of E; can be 
determined uniquely by this vertex through the complementarity interchanges of the 
coordinate values. €; N €; = @, if i # j. The procedure is called “cube-splitting” of 37,41. 

(2) The “cube-splitting” of 27,,, keeps the monotonicity property in the following way: let F be a 
monotone binary function defined on £/1,,, then either Np N €; is empty, or it satisfies the 
binary monotonicity property, 1.e., for arbitrary vertex a of Nz N &;, all vertices of E; greater 
than a, also belong to Ne N €; (fori = 1,--:,n). 

By integrating (1) and (2), a novel monotone recognition method has been proposed in [6]. 


2.1 Definitions/descriptions 


For each vertex V; = (Vj,,°"+, V;,,) of H we compose the following set: 
Ey, = (Gi An) © Fins ila; € {Vi,,m — v;,} for all j, 1 <j < n}, 
and call Ey, the vertical equivalence class of V;. 
Ey, contains a unique vertex from Hf, - this is the vertex with all coordinates > m/2; and 
contains a unique vertex from H, - this is the vertex with all coordinates < m/2. The remaining 
vertices of Ey, can be obtained by component value inversions (with respect to m). Ey, N Ey, = 
® for different vertices V; and V; of H. In this manner 5”, can be split into |A | disjoint sets 


/equivalence classes/ uniquely defined by the elements of A (or H). 
The number of elements of Ey, varies between 2° and 2” depending on the number of 


components of V; differing from m/2. Indeed, if k denotes the number of components of V; 
differing from m/2, ie. k = |{vj,|v;, # (m — v;,)}|, then |Ey,| = 2*. Notice that k = n always 


for odd m. For example, Fig. 2 demonstrates E(3 4,3), €(2,3,4) and E(4,2,2) in Be. 










‘2 
Y 
Me 
i 
K/ 


V7 

x 

Ae 
e3 

“ 


i ‘ 
& 
‘ 


\) 

‘a 

x \ 
e 
ye 

/ 
Ne 
KY 
\ 
t 







aR 
KN 
Ag 
S 


( 
\) 
.) 
A\ 
NS 








fe 
aC 
vi 
Ng 
VA 
Of 
KEK, 
Me, 
\8 
(\? 
Xx) 
VB 
4 









: 
; 
i 
LY 
) 
v 









(4 
9 
Y 

f 


y 
NE 
Mie 
e 
@ 
Vis 
\ 
) 
Xd 
V5 
LX 4 
Xe 
iV 
\ 
d 


)sx 






Z) 
7 
.) 
\ 
i) 

Ag 

/\ 
ay 
Lae 
‘a 
le 
& 

@) 





y 


(2) 
Xo 
\ 
0): 
ue 
NV 
LE 
AX 
XY 
Wi 
hy 
OY 
hE 
hy 
\ 
LX 


V 





r 
5 
Ag 
Nx 
(a 
Ne 
NB 
NB 










VJ 
fy 
War 
1G 
Vy 
ws 
ie 
®  \e 
V7 
\ 





Fig. 2. Examples of cubes in a cube split of 22. 
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For every vertex (@1,°++,@,) of the equivalence class Ey,, we distinguish its sub-list of all 
coordinates accepting a value differing from m/2, let this be the list: (a,,,---,@s,). This list 
exactly fits the list of all coordinates of V; that are different from m/2. The reminder part of 
coordinates accepts the only value m/2 over the V;, as well as over the whole set of vertices of 
Ey,. Now, identifying (a,,,-*-,ds,) with the binary sequence 6 = (f4,---,6,) of length k such 
that 6; = 1 if and only if As; > m/2, - we map (@j,°*:,@,) into the vertex B = (f1,°::, Px) of 
the k-dimensional binary cube E*. In this manner, we obtain a 1-1 mapping M: ey, SE k The 
vertex of Ey, with all coordinates < m/2 is mapped into the vertex (0,---,0) of E K (on the 0-th 
layer); the vertex of Ey, with all coordinates > m/2 is mapped into the vertex (1,---,1) of E e 
(on the n-th layer); and, in general, all vertices of €y,, which have J coordinates > m/2 
(consequently, m — 1 coordinates < m/2), are mapped into the vertices of I-th layer of E*. 

Hereafter, all structures (vertices, chains, cubes, functions, etc.) in 27 ,, will be referred 
to as origin; and all structures (vertices, chains, cubes, functions, etc.) in binary cubes will be 
referred to as induced. 

For example, the induced binary cubes for €(3,4,3), E(2,3,4) and E4,2,2), are given in Figure 
3, (a), (b), and (c), correspondingly. 





(b) (c) 


Fig. 3. 


2.2 The Algorithmic Framework 


Let F be a monotone binary function (which should be recognized with the help of an oracle), 
defined on =" ,,, and let Nr denote its set of units of function F. 


Algorithm 1 
In a theoretical level description, the algorithm implements the following steps: 

1. Apply the cube-splitting on 57,41 and let €),---,.€)q, be the equivalence classes of the upper 
homogeneous elements. 
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2. Compose the corresponding induced binary cube EF; for every €;. 

3. Apply the Hansel’s algorithm to recognize the induced Boolean function f;, defined on E; as 
follows: for every B € E;, f,(8) = 1 if and only if F(b) = 1, where b is the origin of # in 
en, (fori =1,--, |A|). 

4. Transfer the recognition results into the £7), 4. 


3. Implementation 


3.1 Implementation details of Step 1 and Step2 


We consider the lexicographic ordering of the vertices of |A F where the smallest numerical 
values of coordinates are coming first. Thus, the smallest vertex of H in this ordering is 
ce, a) if m is odd, and it is the vertex(>, =, 5), if m is even; the greatest 
vertex is (m,m,:,m). Henceforth, we will assume that A = {V,,V>,-,Vq} is the 
lexicographically ordered set of upper homogeneous elements. 

The cube splitting of 27,4, assumes that we compose for every vertex V; = (Vj,,°-, ¥j,,) of A 
its vertical equivalence class Ey,, and the corresponding induced binary cube to this. 

But at this point we do not need to compose and keep (and further to map to the binary cube) 
the whole set €y,; instead, with every vertex V; = (¥;,,°+:, Vi,,) of Hf we will keep the following 
parameters: 

- the number Ty, of coordinates of V; differing from m/2, - this will determine the size of the 
induced binary cube E;. When V;, that is the issue, is evident, we will just use the notion Tt for 
this, 

- the positions of coordinates differing from m/2, we denote it by the vector V;,, and 

- the values of coordinates differing from m/2, we denote as the vector Viz. 

T, Viz, and Viy will allow the easy reverse mapping, RM : E* — €y,, 1.e., will allow to recover 
Ey,- 

For example, with the vertex (2,3,4) of 22 we keep: 

- numerical value 2, - this is the number of its coordinates differing from m/2, 

- indexes 2,3 - these are the coordinate indexes, where the values are differing from m/2, 

- and 3,4 are the values at the coordinates 2 and 3. 

E(2,3,4) is mapped into the 2-dimensional binary cube E? according to 2"¢ and 3" coordinates of 
(2,3,4), and the accompanying vectors are V;z = (2,3) and Viz = (3,4). 

The reverse mapping is as follows: given the pair (2,3,4) and E? or alternatively, Viz, Viy and 
E?. Consider an arbitrary vertex of E?, for example, let 8 = (1,0), then it follows that the origin 
of B in E(2,3,4) is € = (2,3,0), because: 

- the first component, missing at V;,, must be m/2, that is, a, = 2, 

- the second component should not be inverted in accord to B = (1,0), and, thus, a, = 3, 

- the third component should be inverted, and thus, az = 4-4 = 0. 

In general, let (8,,--:,8;,) be an arbitrary vertex of the k-dimensional binary cube E*, and let 
(V1, °*+, Un) be the upper homogeneous vector of the origin of E*, and suppose that s,,---, Sz are 
its coordinates differing from m/2. Then, the origin of (B,,°*:, By) is (a4, °**, An), where: 
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Vs; if B; =1 
m= vs if fb; =0 
a; = m/2 for i # $y,°*+, Sx. 


as, = 


for j =1,-+,k (1) 


3.2. Implementation Details of Step 3 and Step 4 


In this part, Algorithm 1 recognizes monotone Boolean functions in the |A | number of binary 
cubes of different sizes (some of the functions might be identically 0, but we do not know this 
fact beforehand), by applying the Hansel’s algorithm. 

Also at this step, we will not deal with the binary cubes themselves, and we will be using 
the chain algebras, and therefore, we have to map (by the reverse mapping) all induced structures 
(vertices, functions, chains, etc.) into their origins in £7, ,,. 

For example, consider some monotone Boolean function f; on a k-dimensional binary 
cube. If we obtain the value f;(8) on some vertex # (in the process of the Hansel’s algorithm), 
and know its origin upper homogeneous vertex (V,°*:,V,) and also its coordinates differing 
from m/2 (let they be s,-+-, s,), then we can set that F(@) = f,(@), where the coordinate values 
of @ are defined in accord to (1). 

Similarly, we can map chains from the induced binary cubes into the 27,4. 

For example, the maximum length chain <(000), (100), (110), (111)> in the cube in 
Figure 2 (a) (which is induced to €(343)), 1s mapped to the origin chain <(101), (301), (341), 
(343)> in 2744. 

Upon receipt of the oracle’s response for a given vertex of the binary cube - the response 
is mapped into the origin vertex of 2/7, ,,. Certainly, the response value could also be extended 
by the monotonicity property to the other relevant vertices of £7,,,. But in this research we 
prefer and emphasize the opportunity of the parallel implementation of the recognition 
algorithms in all the induced binary cubes, and so we keep them as separate nonintersecting 
processes. 


4. Small Number of Units 


Note that Algorithm 1 is worth applying when a large number of unit vertices of the monotone 

function appear in the upper homogeneous area H. 

However, in some cases, mostly coming from applications, the function to be recognized has 

a small number of unit vertices in the upper area A, or its complement has a small number of 

zero vertices in the lower homogeneous area H. In the latter case we can recognize the 

complement of the monotone function. For this reason, the following points will be taken into 
account: 

- &y, can be defined for each vertex of H (obviously we will obtain the same set). In the 
example given in Figure 2 the highlighted sets of H will demonstrate E(1,0,1)> €(2,1,0) and 
E(o,2,2), a8 well. 

- The mapping M:é), > E k will be defined as follows: for every vertex (a1,°**, dn) of Ey, 
let (ds,,°**,@s,) denote its subsequence with all coordinates differing from m/2. Identify 
(as,,°+, @s,) with the binary sequence 6 = (f,,-+-,B,) of length k such that 6; = 1 if and 
only if As, < m/2. The vertex of €y, with all coordinate values < m/2 is mapped into the 
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vertex (1,---,1) of E*, - this is on the n-th layer; the vertex of Ey, with all coordinate values 
< m/2 is mapped into the vertex (0,---,0) of E*, - this is on the 0-th layer; and, in general, 
all vertices of Ey, which have | coordinates < m/2 (consequently, m —1 coordinates > 
m/2) are mapped into the vertices of l-th layer of E*. 

- The reverse mapping: RM: E; — &y, will be implemented as follows: let (Aj,---,B,) be an 
arbitrary vertex of the k-dimensional binary cube E*, and let (v,,--+,U,) be the lower 


homogeneous vector of the origin of E*, where s,,+++,S, are coordinates differing from m/ 
2. 


Then the origin of (61,°:+, Bx) iS (@y,°**, An), where: 


for j =1,--,k (2) 


J 


i Vs; if B; = 0 
a m—v, if Bj =1 
a; = m/2 for i # $4,°*+, Sx. 


4.1 Constraints 


Consider the application, which is the generalized model of the known association rule mining - 
in case where in addition to the presence or absence of elements in itemsets, the number of their 
repetitions is also included. The details are given in [7]. Here we highlight the following 
constraints/restrictions that may appear with this problem. In terminology of supermarket basket 
analysis, here we distinguish two postulations: some item exists in the current basket, and 
second, which is the actual number of that item in the basket. 

Let a; be the repetition number of the i-th element, for i = 1,---,n. 

(1) the classic case is the (0,1) vector of item indicators in baskets, basket inventory. 

(2) ay +++ +@, <7, - the summary number of elements’ repetitions (the basket volume) is 
restricted by r, 

(3) a; < 7;, -the repetition number of each aq; is restricted by 7;, the item purchase restriction. 


In these cases, the problem deals with the recognition of monotone functions, where: 
(2) the zeros of the function appear in lower layers of the multivalued grid, 
(3) the zeros of the function appear in some homogeneous bottom area of the grid. 


In both cases it is more efficient to use the second algorithmic scheme of [6]. 
The idea is as follows. Let F be a monotone function defined on 2”, ,. First we note that Nr N H 
satisfies the monotonicity property, i.e., for arbitrary a,b of H, ifa > b then F(a) => F(b). 


Algorithm 2 
1. Firstly identify the part of the monotone function belonging to H by one of the known 
resources of identification of monotone functions, and thus, reduce the size of the multi- 
valued grid. As a result we obtain N; NH. 
2. Apply the cube-splitting according to Nz NH, that considers the vertical equivalence 
classes €j,-++,.E)y,.qq| only for the vertices of Np N A. 
3. Implement 2-4 Steps of Algorithm 1. 


The algorithm can easily be adjusted for the identification of the complement of F in £7) ,,. 
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4.2 Resources 


To implement Step 1 of Algorithm 2 we have the following resources of identification of 
monotone functions in 7. 

a) The first resource is the known algorithm by V. Alexeyev [2]. We notice that applying 
the algorithm to identification of monotone functions on A (instead of 2" 4,4) is a significant 


reduction of the work. The complexity of the algorithm Up on Af will be: 
T (Up) < |My| + |logaC)| - INil, 
where M, and N, are the middle layers of H, defined accordingly. 


b) — The second resource that we may use, is - applying the cube-splitting itself for identifying 


the function on H. We define the upper homogeneous area A of H,; this is the set of all elements 
of £4, with all coordinate values > 3m/4; then 


1. Apply the cube-splitting on A and find €,,-++,.€ A> equivalence classes according to the 


elements of A. 

2. Compose the corresponding induced binary cube E; for every €;. 

3. Apply the Hansel’s algorithm to recognize the induced Boolean function f;, defined on 
E; as follows: for every B € E;, f,(8) =1 if and only if F(b) = 1, where b is the 
origin of B in A (fori = 1,--, |A)). 

4. Transfer the recognition results to A. 

At this point we find Ng N A, and then continue with: 

5. Apply the cube-splitting according to Ne NA, that is find the vertical equivalence 
classes E,,++,.E|y,,nq| Only for the vertices of Np N A. 

6. Implement 2-4 Steps of Algorithm 1. 


c) The cube-splitting can be applied recursively. 
The following resource is also worth mentioning that can be used in all cases/algorithms. This is 
the growing technique [8] in monotone Boolean function recognition and chain computation 
algebra [9]. 


5. Conclusion 


The problem of query based algorithmic recognition of monotone binary functions defined in 
multi-valued grids is known as a hard problem. It was investigated by V. Korobkov, V. 
Alekseev, A.Serjantov, and others. It is known a chain-split type algorithm developed by V. 
Alekseev, where the complexity estimates were given in terms of sizes of the middle layers of 
the grid. A novel method of identification of monotone binary functions based on the partitioning 
of the grid into discrete structures isomorphic to binary cubes was proposed in our recent work, 
where a theoretical level description of two algorithmic schemes was introduced. This paper 
provides the implementation details of the algorithms, as well as focuses on the recognition of 
monotone binary functions with a small number of units /or zeros/ distributed in homogeneous 
areas of the grid. 
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AHHOTanHA 


Hospiii 1o7X0, MOHOTOHHOTO paciio3HaBaHHa Ha OCHOBe pa30HeHHA MHOTO3Ha4HON pellieTKU 
Ha JIMCKpeTHBIe CTpyKTypbI, H30MOpHbIe OuHAapHbIM KyOaM (MeTOX “KyOuYecKoro pa30neHHA’’) 
Iipeq1O%KeH B CepHH TocueqHUX padoT, rye Ha TeopeTH4ecKOM ypoBHe JaHO OMMcaHHe JBYyX 
QIFOpHTMoOB /ayIropHTMH4ecKHxX cxeM/ pellieHua 3afa4uu. B jaHHOol cTaTbe TpHBOANTCA 
ToApoOuHoe onMcaHve JeTaet peawu3alHu ITUX aITOPHTMOB, a TaKKe paccMaTpHBaeTcaA cydali 
pacrio3HaBaHHA MOHOTOHHOM PyHKUHH C HeOOJbIIMM YHCIIOM CMH, YTO CBA3AHO C PAJIOM 
IIpakTH4eCKHX MpHIO*KeHHH. 

KunoueBbie cJ10Ba: paciio3HaBaHie MOHOTOHHOM dyHKUWH, MHOro3Ha4Had pelllerkKa, 
KyOu4eckoe pa30ueHne. 


